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ABSTRACT 

We present the critical conditions for hot trans-fast magnetohydrodynamical (MHD) flows in a 
stationary and axisymmetric black-hole magnetosphere. To accrete onto the black hole, the MHD 
flow injected from a plasma source with low velocity must pass through the fast magnetosonic 
point after passing through the "inner" or "outer" Alfven point. We find that a trans-fast MHD 
accretion solution related to the inner Alfven point is invalid when the hydrodynamical effects on 
the MHD flow dominate at the magnetosonic point, while the other accretion solution related to 
the outer Alfven point is invalid when the total angular momentum of the MHD flow is seriously 
large. When both regimes of the accretion solutions are valid in the black hole magnetosphere, 
we can expect the transition between the two regimes. The variety of these solutions would be 
important in many highly energetic astrophysical situations. 

Subject headings: accretion — black hole physics — MHD — relativity 

1. Introduction 

In order to explain the activity of active galactic nuclei (AGNs) and compact X-ray sources, we consider 
a black hole magnetosphere in the center of these objects. The magnetosphere is composed of a central black 
hole with surrounding plasmas and a large scale magnetic fleld. The magnetic held is originated from an 
accretion disk rotating around the black hole. The electrodynamics of the black hole magnetosphere has 
been discussed by many authors; force-free magnetospheres were discussed in Thorne, Price, & Macdonald 
(1986) and more general magnetospheres in Punsly (2001). 

In the black hole magnetosphere, because of the strong gravity of the black hole and the rapid rotation 
of the magnetic field, both an ingoing plasma flow (accretion) and an accelerated outgoing plasma (wind /jet) 
should be created. The plasma would be provided from the disk surface and its corona. When the plasma 
density in the magnetosphere is somewhat large, the plasma inertia effects should be important. In this 
case, the plasma would be nearly neutral and should be treated by the ideal magnetohydrodynamic (MHD) 
approximation (Phinney 1983), so the plasma streams along a magnetic field line, where the magnetic field 
line could extend from the disk surface to the event horizon or a far distant region (Nitta, Takahashi, 
& Tomimatsu (1991); see also Tomimatsu & Takahashi (2001)). The outgoing fiow effectively carries the 
angular momentum from the plasma source, and then the accretion would continue to be stationary, releasing 
its gravitational energy. The magnetic field lines connecting the black hole with the disk, which are mainly 
generated by the disk current, may not connect directly to the distant region, but via the disk's interior the 
energy and angular momentum of the black hole can be carried to the distant region; the energy and angular 
momentum transport inside the disk is not discussed here. 
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If the plasma density is sufBciently low and the magnetosphere is magnetically dominated, one can 
expect the pair-production region along open magnetic flux tubes, which connect the black hole to far 
distant regions directly (Beskin 1997; Punsly 2001). Further, we also expect the Blandford-Znajek (1977) 
process, which suggests the extraction of energy and angular momentum from the spinning black hole. 

In this paper, we assume a stationary and axisymmetric magnetosphere, and consider ideal MHD flows 
along a magnetic field line. The initial velocity can be at most less than the slow magnetosonic wave speed. 
To accrete onto the black hole, the ejected inflows from the the plasma source must pass through the slow 
magnetosonic point (S), the Alfven point (A) and the fast magnetosonic point (F) in this order, as it is 
well known. At these points. A, F and S, the poloidal velocity equals one of the Alfven wave and fast and 
slow magnetosonic wave speeds, respectively. In the case of accretion onto a star, because the accreting 
plasma is stopped at the stellar surface, a shock front would be formed somewhere on the way to the 
stellar surface and the accretion becomes sub-fast magnetosonic. However, for accretion onto a black hole, 
the flow must be super- fast magnetosonic at the event horizon (H). If not so, the fast magnetosonic wave 
can extract information from the interior of the black hole to the exterior; this fact obviously contradicts 
with the definition of the event horizon. In fact, an ideal MHD accretion solution which keeps sub- fast 
magnetosonic has zero poloidal velocity at the event horizon and the density of the plasma diverges; the 
solution is unphysical. 

Because the magnetic field lines would rigidly rotate under the ideal MHD assumption, there are two 

light surfaces (L) in the black hole magnetosphere (Znajek 1977; Takahashi et al. 1990, hereafter Paper I). 
The plasma source must be located between these two surfaces. Further, one or two Alfven surfaces lie 
between the two light surfaces, and for accretion there must be a fast-magnetosonic surface between the 
Alfven surface and the event horizon (see Paper I). Here, we should note that the physical mechanism to 
determine the angular velocity of the field lines is controversial. A time-dependent determination of it has 
been discussed by Punsly (2001); the torsional Alfven wave originated from the plasma source and propagated 
up and down the magnetic flux tube forces to minimize the magnetic stresses in the system. 

The conditions on the flows at the magnetosonic points and the Alfven point restrict the five physical 
parameters which specify the fiow (sec the following section) if one is to adhere to the ideal MHD assumption 
globally. The fast and slow critical points have X-type (physical) or 0-type (unphysical) topology for the 
solution, while the Alfven point does not specify its topological feature; hereafter, we will call the fiow 
passing through both X-type fast and X-type slow magnetosonic points as the "SAF-solution" . When we 
discuss the global features of a solution, it is very important to know the numbers of these critical points and 
the Alfven points. In a zero- temperature limit (cold limit), the regularity condition for the trans- fast MHD 
flow was discussed by Takahashi (1994). In this case, the Alfven points and the fast magnetosonic points 
only appear in wind and accretion solutions without the slow magnetosonic point, because the velocity of 
the slow magnetosonic wave speed is zero. The rclativistic hot MHD flow equation has been formulated by 
Camenzind (1986a,b, 1987, 1989); see also, Paper 1. In § 2, we summarize the basic equations for the MHD 
flows and the condition at the Alfven point discussed in Paper I. 

Though we consider the general rclativistic plasma flow, the significance of the Alfven point and the 
fast and slow magnetosonic point conditions is similar to that of a Newtonian wind model by Weber & Davis 
(1967) and a special rclativistic wind model by Kennel, Fujimura & Okamoto (1983). Kennel, Fujimura 
& Okamoto (1983) classified the outgoing trans-Alfvenic MHD wind solutions into a "critical" {E = Ep) 
solution, "sub-critical" {E < Ep) solutions and "super-critical" {E > £^f) solutions, where E is the conserved 
energy of the wind and Ep is the energy for the trans-fast MHD wind. To reach distant regions, the critical 
(trans-fast MHD wind) solution and the super-critical (sub-fast MHD winds) solutions are physical, and the 
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sub-critical solutions are unphysical beyond the turnaround point. Turning now to the ingoing plasma flow, 
the topology of the black hole accretion solution space also has a similar structure; that is, (i) a trans- Alfven 
MHD ingoing flow with E = E-p, which means a trans- fast MHD ingoing flow discussed in this paper (critical), 
(ii) trans-Alfven MHD ingoing flows with E < Ep (sub-critical) and (iii) trans-Alfven MHD ingoing flows 
with E > E-p (super-critical) . Under the ideal MHD approximation, the physical solution is only the critical 
solution (i); the super-critical solutions (iii) are unphysical for the reason mentioned above. In addition to 
these, for accretion onto a black hole, we must consider (iv) sub- Alfven (or sub-slow MHD) ingoing flows, 
although they do not pass through the Alfven point A. The breakdown of ideal MHD approximation between 
the horizon and the inner light surface is indicated by Punsly (2001), and then non- ideal MHD solutions 
classified into (ii), (iii) and (iv) would be realized as accretion solutions onto a black hole (sec § 5). 

The main purpose of this paper is to examine the thermal effects on an ideal MHD plasma streaming 
in a black hole magnetosphere (see Fig. 1) by studying the critical conditions at those magnetosonic points. 
Now, the slow magnetosonic point appears on the MHD flow solutions. The details of critical conditions 
at the fast and slow magnetosonic points arc discussed in § 3. We derive the critical conditions at the fast 
and slow magnetosonic points, which are denoted in terms of the location of the fast and slow magnetosonic 
points, the sound velocity at the fast and slow magnetosonic points and the locations of the Alfven and light 
surfaces. In § 4, we clarify the thermal effects on the MHD flows, and discuss its dependence on the rotation 
of the black hole magnetosphere and the divergence of the cross-section of a magnetic flux-tube along the 
field line. Then, we can find two kinds of trans-fast MHD fiow solutions for both inflows and outflows: 
"hydro-like" MHD flow and "magneto-like" MHD flow. The main difference between the two solutions is 
the behavior of the magnetization parameter, which is the ratio of the fluid and electromagnetic parts of the 
total energy of the flow. The hydro-like MHD flow solution is a somewhat hydrodynamical solution and, in 
the weak magnetic field limit, this trans- fast magnetosonic flow solution becomes a trans-sonic flow solution 
discussed by Abramowicz (1981) and Lu (1986) in the hydrodynamical case. We also unify hydrodynamic 
flows with hot MHD flows in a common formalism. The hydro-like MHD flow solution disappears for a 
magnetically-dominated magnetosphere. On the contrary, the magneto-like MHD flow solution results in 
the magnetically-dominated flow, although that disappears for hotter plasma cases. In § 5, we summarize 
our results. 



Wc present basic equations of a stationary and axisymmetric ideal MHD flow. The flow streams along 
a magnetic field line in the black hole magnetosphere, and accretes onto the black hole or blows away to 
a far distant region. To determine the configuration of magnetic field lines and the velocity of MHD flows 
streaming along each magnetic field line, we must solve self-consistently what is called the Bernoulli equation 
along magnetic field lines and the magnetic force-balance equation. These equations are derived from the 
equation of motion for relativistic MHD plasma 



the conservation law for particle number [nu^).^ = 0, the ideal MHD condition u^F^^j, = and Maxwell's 
equations. Here, p, P and n are the total energy density, the pressure of the plasma and the proper particle 
number density. The electromagnetic field tensor F'^" satisfies Maxwell equations and is the four- velocity 



2. Basic Equations and Trans- Alfvenic MHD Flows 




(1) 
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of the plasma. The background metric is written by the Boyer-Lindquist coordinates with c = G = 1, 
, o 2mr\ ,0 ^amrsiv?6 , 

ds^ = ( 1 — j dt^ + dtd(t) (2) 

+ af ^ = sin^ Odcjf - —dr^ - UdO^ , 



where A = — 2mr + a^, S = + cos^ ^, and m and a denote the mass and angular momentum per 
unit mass of the black hole, respectively. 

The flow with the above assumptions streams along a magnetic field line, which is expressed by a 
magnetic stream function = '^{r,6) with — constant; ^ is basically the toroidal component of the 
vector potential. Stationarity, axisymmetry and ideal MHD condition require the existence of five constants 
of motion (e.g., Bekenstein & Oron 1978; Camenzind 1986a). These conserved quantities are the total energy 
E{'^), the total angular momentum L{^), the angular velocity of the field line f2i?(^') and the particle flux 
through a flux tube which are given by 

L = -^U4, - -^B^ , (4) 

= -TT- = --FT- , (5) 

V = -5^ , (6) 

-Dp 

where B^ = (A/S) sin 6* Fgr is the toroidal component of the magnetic field. Bp is the poloidal component 
of the magnetic field seen by a lab-frame observer 

Bp = [g'-'-idr^f+g'^ide^iff] , (7) 

Pw 

and = gf^ — gttg4>4>- The poloidal component Up of the velocity is defined by Up = —uau^ {A = r, 0), 
where we set Up > {u"^ < 0) for ingoing flows. For the polytropic equation of state with adiabatic index F, 
the relativistic speciflc enthalpy ^ is written as (see Camenzind 1987), 



/U = rrip 

where 



UpBp^J 



(8) 



'''"^-F-lninimp ' 

and nip is the rest mass of the particle. The quantities labeled by "inj" are specified at a point injecting 

plasma as a plasma source. The boundary conditions should be given by a plasma source model (e.g., the 
accretion disk/corona model, pair-creation model, and so on). Thus, we require the specification of the 
additional fifth constant of motion, hiny, in a cold limit, we obtain hinj — > and /x nip (= ^c)- 

By using the conserved quantities, the equation of motion projected onto the direction of a poloidal 
magnetic field, which is called the poloidal equation (and is often referred to as the relativistic Bernoulli 
equation), can be expressed by (e.g.. Paper I) 

{l+ul) = {E/t,f[{a-2M')f-k] , (10) 
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where 



a = gtt + '^9t4>^F + 9<t,<t>^]r > (11) 
k = {g^^ + 2gt^L + gttL^)/pl , (12) 
f ^ _ (gt0 + 9<t"t'^F) + {gtt + 9t^^F)L 

and L = L/E. The relativistic Alfven Mach- number M is defined by 

^2 ^ ^T^l^nul ^ 4wfir]Up _ ^^^^ 
Bp Bp 

Note that a~^/^ is the "gravitational Lorentz factor" of the plasma rotating with the angular velocity flp in 

the Kerr geometry, whose definition includes both the effects of the gravitational red-shift and the relativistic 
bulk motion in the toroidal direction. The locations of the Alfven points (rA, 0a) along a magnetic field line, 
where 6 = 6{r; are defined by = a. 

The relativistic force-balance equation, which is the equation of motion projected perpendicular to the 
magnetic surfaces, was derived by Nitta, Takahashi, & Tomimatsu (1991); see also Beskin (1997). We should 

solve both the poloidal equation and the force-balance equation, but it is too difficult to solve these equations 
self-consistently. So, we only discuss the poloidal equation on a given magnetic field line. When the poloidal 
field geometry ^ = ^(r, 9) is known and the various conserved quantities are specified at the injection point, 
the Mach-number (14) used together with equation(8) in the poloidal equation (10) determines a complicated 
equation for Up as a function of r. It seems that, to obtain the poloidal velocity, we must solve a polynomial 
of degree 4A'' + 2 in z = Up^^ for the polytropic index N = l/(r — 1) (Camenzind 1987). To study the 
behavior of accretion and wind/jet solutions, however, we can reduce the poloidal equation to 

where TZ = e^a — 2e^M^ — kM^ and e = 1 — QpL, and we can plot a contour map of {E/ easily on the 
r-Up plane under the given flow parameters Q.f, L, t] and /linj- The regions with 7^ < 0, where E"^ < and 
no physical flow solution exists, give the forbidden regions of the flow solution on the r-Up plane; the shape 
of the forbidden regions are classified hy ^Ip and L (see Paper I for the details of the forbidden regions). 
Note that, near the inner and outer light surfaces (r = and r = r£"*) defined by a = 0, the toroidal 
velocity of the plasma approaches the light speed. However, if the ideal MHD condition breaks there because 
of the inertia effects, the plasma would enter the forbidden region _R < by crossing the magnetic field lines; 
the non-ideal MHD flows are no longer forbidden in the R < regions. The breakdown of ideal MHD for 
accretion near the inner light surface is discussed in the last section. 

At the Alfven point (A), where r = rx and = a, it seems that the function / diverges. However, 
to obtain a physical accretion solution passing through this point smoothly, we must require the following 
condition 

l^tM^^Z^^ (16) 

where uj = —914,/ 9(b<t> is the angular velocity of the zero angular momentum observer with respect to distant 
observers and the subscript "A" means the quantities at the Alfven point. Thus, the ratio of the total 
angular momentum to the total energy of the flow is determined by the location of the Alfven point and 
Vtp. When Q < Q,p < uj^ (classified as "type H" in Paper I), where wh is the angular velocity of the 
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black hole, one Alfven point A™ or A°"* appears in the r-Up plane. When loh < < fimax ("type I") 
or ^^min < ftp < ("type HI"), two Alfven points A"* and A™/°"* appear, where Omin and f2max are the 
minimum and maximum angular frequency for the light surfaces to exist in the magnetosphere. If two Alfven 
points appear in the r-Up plane, the physical flow solution passes through one of them. The Alfven point 
labeled by "m" or "out" corresponds to inside or outside the separation point r = Tgp (SP) which is defined 
by a' = 0. At r = Tsp, the gravitational force and the centrifugal force on plasma are balanced if the poloidal 
velocity of the plasma is zero, hence the term separation point. Although the Alfven point is classified 
into three types, one can interpret that the MHD flow of type I, II or III accretes onto a slow-rotating, 
rapid-rotating or counter-rotating black hole, respectively, when the value oiflp is specifled. 

From equations (10) and (16), we can express the total energy E and the total angular momentum L 
as functions of the Alfven radius and the injection point as follows: 

^ ^ {9tt+ 9t<t,^F)A c 

Q-A 

L = i-94>4>)^{^F-^A)^ ^ .^g. 

"A 

where £ = E — flpL > can be specified at the plasma injection point and the Alfven point. The condition 
of a negative energy MHD accretion flow, {gu + gt<p^F)A < 0, is unchanged from the cold hmit case (see 
Paper I). Thermal eflFects on the hot plasma flow are only included in £, and modify the amplitude of the 
ingoing energy and angular momentum. 



3. Critical Conditions at Magnetosonic Points 

For a physical solution for accretion onto a black hole, we also require that the critical conditions at 
both fast and slow magnetosonic points must be satisfied. In this section, we discuss restrictions on the 
remaining field-aligned parameters; we see that the locations of the fast and slow magnetosonic points give 
the total energy and the particle number flux along the magnetic field lines. 



The differential form of the poloidal equation (10) is written by (see Paper I) 

V 



(Inup)' = - , (19) 



where 



^^ = [f) [ [^(^' - ")Csw + MU^] {InBpY + ^(l + C^J [M\M^ - a)k' - Qa'] | , (20) 

V = {M^-af[{C^^-ul){M^-a) + {l+ul)MWTZ-'] (21) 

with A'^ = e'^ + ak = f{M'^ - af and Q = ae^- Se^M^ - 2kM'^. The prime (• • •)' denotes [{d0^)dr - 
{dr^)do]/{^/^Bp) which is a derivative along a stream line. The relativistic sound velocity Osw is given by 

and the sound four- velocity is given by C^^ = a^^/ (1 — a^^). 
The denominator (21) can be reduced to the form 

Voc{ul- uiw)^ - w|m) (wp - Usm) . (23) 
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where the relativistic Alfven wave speed uaw, the fast magnetosonic wave speed ufm and the slow magne- 
tosonic wave speed usm are defined by 

"^w(^;*) - i^"- (24) 

4m (r-;*) = l(^Z + sjz^ - AC^^ul^^ , (25) 
uluir;^) = l(^Z-^Z^-4Ci^ui^^ , (26) 



with 



= -AnriEp^f . (28) 

When Up = = ^fm O'" = '"smj the denominator becomes zero. Therefore, at these singular points, 

we must require M = Q to obtain physical accretion solutions which pass through these points smoothly. 
The location of = u\ [= Up^{r\]'^)] is the Alfven point discussed in the previous section. Similarly, 
the locations of = Up [= Wfm(''f; ^)] and Up = Ug [= 'Usm(''s; *)] correspond to the fast magnetosonic 
point r = rp and the slow magnetosonic point r = rs, respectively. We should mention that, to calculate the 
Alfven velocity, we need to solve a polynomial of high degree, while in the cold limit it is simply obtained 
as u^L*'^ = (a-Bp)A/(47r^c??); the Alfven velocity of a hot MHD flow is always smaller than u'^^'^. 

Figure 2 shows schematic pictures of the V = curves in the r-u^ plane. [The definition of the poloidal 
velocity Up includes the gravitational-redshift factor, so that Up diverges at the event horizon for a physical 
accretion solution. Hereafter we use the r-u^ plane when we discuss the behavior of accretion solutions; 
under a given magnetic field with ^'(r, 0) = constant, we can also calculate 9 = 9{r; \E') and = u^{r; *). 
] Corresponding to the two modes of magnetosonic wave speeds and the Alfven wave speed, we see three 
branches of D = curves, which correspond to = — '"sm ^'^^ — ^Aw curves. The = u\yf 

curve is always located inside the forbidden region (the shaded regions) except for the Alfven point, so these 
forbidden regions separate the r-w'' plane into one or two super- Alfvenic region(s) and one sub-Alfvenic region. 
If there is no area of k{r; L) > 0, we see one super- Alfvenic and one sub-Alfvenic regions with forbidden 
regions classified as "type A" in Paper I. The total angular momentum has a value of Z_ < L < i+, where 
L = -/j+(> 0) and L = L^{< 0) are the minimum and maximum values of the function L = L{rp^] ^f)- On 
the other hand, if an area with fc > exists between two light surfaces, we see two super- Alfvenic regions 
separated by "type B" forbidden regions, and obtain larger total angular momentum MHD flows (L > -L+ 
or \L\ > Thus, we can classify the forbidden regions by Q,f and L as type I, H or HI and type A or B, 

independently; hereafter, we will denote the type of forbidden regions as, for example, type lA. 

The Up = uIj^.j and = Ugj^.j curves are located in the supcr-Alfvenic region and the sub-Alfvenic 
region, respectively. Figure 2a shows a case for strong magnetic fields satisfying (Cg^)A < u\, and Figure 2b 
shows a case for weak magnetic fields satisfying {C^^)a > u\. In Figure 2a the Up = u^y^ curve connects 
to the Alfven point marked by "A" (r = rA,Up = u\), while in Figure 2b the Up = Ugj^ curve connects to 
the Alfven point. At the Alfven radius r = ta the value of the function / — f{up; va) is zero except at the 
Alfven point A. Then, Up = C^^ is one of the solutions of V{up; va) = (marked by "C" in Fig. 2). A similar 
situation at the (outer) Alfven point can be seen in the Newtonian case (see Heyvaerts & Norman (1989)). 

We also plot a typical curve with J\f = in Figure 3. Crossing of the V = curves and Af = curves in 
the super- or sub-Alfvenic region means the fast or slow magnetosonic point, respectively. In the cold limit. 
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between the Alfven point and the event horizon in the super- Alfvcnic region, crossing of the Af = curve 
and T> = curve always exists regardless of the rj value (Paper I). However, in the case of {C^^)a > u\, 
we cannot find any reason for crossing of these lines. In fact, we find a restriction on the hot trans-fast 
MHD accretion by the thermal effects. In this case, the condition A/' = "D = is not achieved between the 
inner Alfven radius and the event horizon; that is, no physical trans-fast MHD accretion solution exists. 
When the thermal effects dominate over the magnetic effects, we can find that the crossing of these lines 
is only available for smaller \r]\, while for larger |?7| it becomes impossible to generate a physical trans- fast 
MHD accretion solution (see below). In the cold limit, because the slow magnetosonic wave speed is zero, 
the X-type slow magnetosonic point is located just on the r-axis (u^ = line), which is just the separation 
point. When thermal effects are effective, we can see X-type slow magnetosonic points with a sub-slow 
magnetosonic region in the sub-Alfvenic region of the r-u'" plane. 

Now, we will discuss the condition for crossings of the Af{r, Up) = and T>{r, Up) = curves. We use the 
indices "F" and "S" to denote the quantities evaluated at the fast and slow magnetosonic points, respectively, 
and use the index "cr" to unite the quantities at these magnetosonic points. In the following equations, to 
discuss MHD flows passing through the fast or slow magnetosonic point, we can replace the subscript "cr" 
by "F" or "S" . From the condition D = at the fast and slow magnetosonic points, the poloidal velocity at 
these critical points v^^ = Up{rcT; ^) is written as 



nC^i^{M^' - a) + M'^A'^ 



(29) 



7e(M2 - a) - M4^2 

and by use of the definition of Mach-number (14), the particle flux through a flux tube t] is determined by 

where the critical Mach-number M^^ = NP^Tct; ^) is obtained as a solution of A/" = 0, which is a cubic 
equation in M^. Thus, we can express as a function of rcr with given parameters Q.f, L, (ag^)cr and 
^f. The total energy of the trans- fast (or trans-slow) MHD flow is also evaluated at the fast (or slow) 
magnetosonic point rp (or rs) by using the poloidal equation (15). Here, we would like to emphasize that 
(Osw)cr is introduced as a parameter for the thermal effects on the trans-fast MHD flows instead of (asw)inj 
(or /linj). The acceptable ranges for (ag^)s and {a^^)-p are restricted by the critical condition (30), which is 
to be realized as trans-slow and trans-fast MHD accretion, respectively. Thus, all boundary conditions can 
be replaced by the Alfven and magnetosonic conditions. The behavior of ry = r?[rcr; ^f, L, {a1^)„] a, 4*] will 
be discussed in the next section. 



4. Thermal Effects on Trans-Fast MHD Flows 

Let us discuss a trans-fast MHD flow in a black hole magnctosphcre. We consider magnetic flux-tubes 
given by dr'^ — and dg'if = Cr~^ sin 6*; that is, the poloidal magnetic field is denoted as i?p(r;\E') = 
(C/VAE) r~^, where S = 5{^) and C = C(5'); S is related to the divergence of a magnetic flux tube. For 
example, 6 = with C(5') = constant means the split monopole magnetic field (Blandford & Znajek 1977). 
Hereafter, we consider a situation where the plasma streams close to the equatorial plane. We expect that 
the qualitative picture is not drastically changed when we leave the equatorial plane and when we consider 
more complicated field geometries. 
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4.1. Restrictions on Trans-Feist MHD Flow Solutions 

We introduce a new parameter xa to specify the Alfven radius by x\ = {r^p — ta)I{tsp — r™'), where 
< a;A < 1- Though, for a given xa, we obtain a value of L, we may find another value for xa giving the 
same L value in the cases with types I and III; that is, we may see two Alfven points inside the separation 
point. Further, wc introduce i) = /ic|'7/C| and = {<'^sw)cr- 

Figures 4a, 5-8 show relations between fj and for various ^^r values under given parameters flp, L, 
S and the spin parameter o; The locations of Tct = rn, rcr = r™^°"*, Vct = r™^°"* and rcr = are marked 
by H, L, A and SP, respectively. Here, we will discuss the magnetosonic points located inside the outer 
light surface, because for accretion the fast magnetosonic point should be located inside the outer Alfven 
point and the slow magnetosonic point should be located between the inner Alfven point and the outer light 
surface. 



4.. 1.1. General Properties 



In the cold limit, there are three branches of ?) = fflr-p) (solid curves with C,cr = 0.0), while there is no 
r) = fj{rs) curve. For hot MHD flows with ^cr = 0.1,0.2 and 0.3, however, both 77 = fi{rF) and 77 = fi{rs) 
curves (dashed curves) exist. When we try to plot a contour map of (E/fic)"^ on the r-u'^ plane as an accretion 
solution (see, e.g.. Fig. 10a), from the 77 vs. r^ diagram we can find the acceptable locations of the fast/slow 
magnetosonic points, which are shown as the crossings of a 77 = constant line and 77 — 77(rcr;Ccr) curves. 
Hereafter, the magnetosonic points located inside the inner Alfven point are labeled as "in", the middle 
magnetosonic points located between two Alfven points are labeled as "mid" and the magnetosonic points 
located outside the outer Alfven point are labeled as "out". In Figures 4a, 5-8, the location of 77(rcr; Ccr) = 



exists between tf = ru and rp 



The location does not depend on the Cf values. For hot MHD flows 



we can see the cases where 77 = 77 (rp) curves and 77 = 77 (rs) curves are connected at rcr = ta (the locations 
marked by "•"). For the fast magnetosonic point located just inside (or outside) the Alfven point, from 
Af =0 we have 

aA + yA\A\'/^+0{A) (31) 



Ml 



where 



In the rp — » ta (>l 



■4 {6e^[a^,^{lnBp)f -2e^a^^^{lnBpyA' + ka'A'y/^ 
■^-\A\ -kalJ\nBpy + {k/2y 

0) limit, from equation (30), the value of 77 is given by 



(32) 



lim fj 



\{Bp/C)a^ 




p{6e2[a2^(lnSp)']2- 


- 2g2a2^(lnSp)'^' + ka'A'} 


_47r(/x//ic). 


A 




k']akA' + 3e2a2^(lnSp)'A;'} 



- 1 



1/2 



(33) 



Then, 77 for hot MHD accretion passing through the inner-fast magnetosonic point has an upper-limit. Thus, 
the thermal effect restricts the acceptable fj values, while in the cold limit the range is given by < r} < 00. 

Figures 4a, 5 8 also show that the 77 = 77(rp ) and 77 = ?7(rp"*) curves shrink down vertically with increasing 
Cf- When (-p > 0.3, the magnetic effect on the plasma can still remain efficient for the flows passing through 



Tp or r™* which gives a trans-fast MHD flow with 77 <C 1. On the other hand, the 77 = 77(rg") and 77 



) 



curves become almost vertical when 77 is at least several times as large as 77* ; that is, the location of the 
slow magnetosonic point is almost independent of fj, and the curves shift toward the inner and outer light 
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surfaces with increasing respectively. We should note that an ideal MHD accretion flow after passing 
through the inner-slow magnetosonic point is impossible, because no Alfven point is located inside this slow 
magnetosonic point; an outflow may be possible after passing through this slow magnetosonic point, the 
inner Alfven point and the middle fast magnetosonic point, in this order. 

In contrast to the inner-magnetosonic points, a trans-fast MHD flow passing through the middle-fast 

magnetosonic point is always available for any fj values. We see that in Figure 4a the location of the middle- 
fast magnetosonic point is insensitive to both 77 and (p; and the location of the middle-slow magnetosonic 
point is insensitive to only f}, while it shifts inward with increasing ^s- We can also see branches for f/ = r}(rp"*) 
for the fast magnetosonic point located outside the outer Alfven point. An outgoing flow from the black hole 
magnctosphcrc should pass through the outer-fast or middle-fast magnetosonic point after passing through 
the outer or inner Alfven point. In this paper, however, we will focus our attention to trans-fast MHD 
accretion onto a black hole, and omit discussions of the outgoing flows from the magnetosphere. 

Figure 4b shows the total energy E/fic of the trans-magnetosonic flows as a function of r^- For a 
fast magnetosonic point giving a smaller fj value, the energy becomes larger. This means that the larger 
energy is shared by smaller numbers of particles. Corresponding to the "maximum" rj value for 77 = fj^r™), a 
"minimum" value of E / ji^ exists for the E = i?(rp ^°"*) curve. For a hotter flow which includes more thermal 
energy, this minimum value becomes larger than the cooler one. A trans-slow MHD flow passing through 
the inner-slow or outer-slow magnetosonic point has also a lower limit for the total energy. Furthermore, for 
the middle-slow magnetosonic point, the acceptable value of E/ is restricted within a very narrow range, 
while it seems that there is no restriction on 7). 

4- 1-2. Rotational Effects of the Magnetic Field Line 

From Figures 5 and 4a, we can see the Clp dependence of the locations of magnetosonic points. Com- 
paring Figure 5a {^f = 0.9f2niax) with Figure 4a {Hp = O.Sfimax), the location of the inner light surface 
moves outward and the locations of the separation point and the outer light surface move inward with in- 
creasing i}p. In both cases, the separation point is located between two Alfven points. For cooler accretion 
of Cf < 0.1, three fast magnetosonic points may be possible between the inner Alfven point and the event 
horizon (see Fig. 5a). The middle one is an 0-type point (imphysical) , while the others are X-type critical 
points (physical). For hotter accretion (^p > 0.1 at least), however, the X-type fast magnetosonic point 
located next to the Alfven point disappears, and it turns to the slow magnetosonic point. 

Next, comparing Figure 5b {Qp = 0.5f2max) with Figure 4a, the location of the inner light surface moves 
inward and the locations of the separation point and the outer light surface move outward with decreasing 

rip- Two Alfven points arc located between the inner light surface and the separation point; and then, both 
the inner and outer Alfven points, which arc labeled by "m" (see § 2), can be related to accretion started near 
the separation point. Concerning the fj = fl{rp"^) branches, which are located between two Alfven points, 
a branch of smaller (p (e.g., Cf = 0.1) has a maximum, while a branch of larger Cf (e.g., Cf = 0.2, 0.3) has 
no upper limit. We see that for smaller Cf the branch for the middle- fast magnetosonic point connects to 
the outer Alfven point (see the Cf = 0-1 curve) and for larger Cf the branch for the outer- fast magnetosonic 
point connects to the outer Alfven point (see Cf = 0.2). 

Figure 6 shows fj as a function of ra with a;A = 0.5. The outer Alfven point is located inside the 
separation point. There is an upper limit to the fj = ?7(rp"^) curve for each of the Cf = 0.1 and 0.2 cases, 
while for Cf = 0.3 the middle- fast magnetosonic point always exists for any fj values (no upper limit). 
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4.1.3. Black Hole 's Spin Effects 

Figures 7a and 7b show the fj vs. Vcr relation with a = 0.5m (a corotating black hole with the mag- 
netosphere) and a = —0.5m (a counter-rotating black hole), respectively. They can be also compared with 
Figure 4a, which is the case with a = 0. In the case with a — 0.5m, two Alfven points arc located inside the 
separation point; and the outer Alfven point is located very close to the separation point. The properties of 
fj = i7(''cr) are similar to the case of Figure 5b. In Figure 7a, it seems that the branches with fj = fj{rp^^) 
have no maximum value, but for cooler flows (Cf 1) there are branches with fj = ^(rp"*) which have a 
maximum value. In the case of a = —0.5m, the separation point is located between two Alfven points, and 
the properties of fj = fi{rci) are essentially similar to the cases for Figures 4a and 5a. Comparing Figure 7b 
with Figure 4, the counter-rotating effect of the black hole generates three (or two) inner-fast magnetosonic 
points, while the effect of corotation is to suppress such multi-inner- fast magnetosonic points generation (see, 
e.g., fj — fi{r^) curves with (-p =0.1 in these figures). 

Figure 8 shows r} as a function of fcr with a = 0.8m, xa = 0.8 and = 0.5f2max- This is a case of 
only one Alfven point and one super- Alfvenic region in the r-u'" plane (type IIA forbidden region). In the 
previous examples of type lA forbidden region (i.e., Figs. 4-8), we have seen inner and outer Alfven points 
and two regimes of trans-fast MHD accretion solutions. In this case, however, there is one type of solution. 
The Alfven point is located inside the separation point, and each branch of fj = fi{rp) has an upper limit for 
fj, except for the cold limit. It seems that for the cooler flows the slow inagnetosoruc point is located near 
the separation point, but for the hotter flows the location shifts outward. However, if the hotter flow has a 
small fj, the slow magnetosonic point could be located close to the separation point. 

4.1.4. Non-conical Effects of the Magnetic Field Geometry 

The effects of magnetic field geometry on the critical condition (30) are shown in Figure 9 for a hot 
MHD flow; The effects of non-conical geometry {S ^ 0) for cold MHD accretion have been discussed by 
Takahashi (1994). The maximum fj value increases with increasing S, for fj = fj{r™) curves. This means that 
the field geometry converging along an ingoing stream line (S > 0) rather than the radial field is available to 
make a higher accretion-rate than ^ < cases. The location of fjirp) = has a weak dependence on 6. The 
location of the middle fast magnetosonic point also has a weak dependence on 6, and the slow magnetosonic 
point remains at a fixed location except for smaller fj cases. 

4.2. Two Regimes of Accretion Solutions 

Here, we will present accretion solutions. To plot an SAF-solution, we need to determine the five field- 
aligned constant quantities: Qp, L, rj, E, hinj. Although these constants should be given as boundary 

conditions at the plasma source (inj), mathematically we can choose the locations of rp, tl, and the 
sound velocity Cf as free parameters to fix the conserved quantities if one restricts their interest to ideal 
MHD solutions. First, for the plasma sources to exist in a black hole magnetosphere, we require that 
f^min < < ^ina,x', then, two light surfaccs r = rf^'-{ilF-.a,'^) and r = r£"*(f7i?; a, 'f) are determined. 
Second, from equation (16), L is determined by the Alfven radius ta, which is located between two light 
surfaces mentioned above. Third, as we have seen in the previous section, by specifying rp and ^f, V and E 
are calculated from equations (30) and (15). Similarly, when rs and Cs are set, rj and E are also calculated. 
Of course, we should require that rj = rj{rF,CF) = '?(?'s,Cs) for the SAF-solution. Finally, we can plot a 



contour map of i?//Xc on the r-u^ plane. The SAF-solution is obtained as the curve with E = Ep = Es, 

where Ecr = E{r„, 

When the forbidden region is type lA (or IIIA), where both inner and outer Alfven points appear on 
the r-u"^ plane, there are two regimes of SAF-MHD accretion solutions (see Takahashi (2000)): 



(i) inj ^ S™^ ^ A'" ^ P" ^ H 

(ii) inj ^ S°"* ^ A™* ^ F™<i ^ H 



The accreting matter for case (i) would be injected near the separation point. In Figure 4b, we see that 
the energy Es/nc is restricted within a very narrow range. Though the sound velocity is not constant along 
the flow, which means Cs 7^ Cf, the possible location of the inner-fast magnetosonic point which must give 
Ep = Eg would be also restricted to a narrow range. Note that if the accreting plasma is intensely heated 
up, there may be no solution of case (i); for example, in Figure 4b, the flow of (s = 0.1 and Cf = 0.2 is 
forbidden, because such plasma heating causes a conflict situation of Ep > Eg. The accreting matter for 
case (ii) would be injected inward from an area between the outer Alfven radius and the outer light surface. 
If the separation point is located outside the outer Alfven point, it is possible that the case (ii) accreting 
matter is injected from near the separation point. Figures 10 and 11 show typical examples of case (i) and 
case (ii), which satisfy the requirement that fj < ?7max(''F )• In Figures 10a and 11a, we also see accretion 
solutions which reach the event horizon with zero radial velocity, but these solutions are unphysical because 
the accreting plasma stops just on the event horizon and its density diverges. 

Now, we introduce Xcm{r) = —ftpBcf,/{4TTri\E\), which means the ratio of the electromagnetic energy 
to the total energy (absolute value); by normalizing with \E\, we will express < — 1 for a negative 

energy {E < 0) inflow. The fluid part of energy per total energy is denoted by E/\E\ — Xem(= -'^fluid), 
where (Xflui(i)H can become negative even if (A'[iuid)inj is positive (Hirotani et al. 1992); that is, the initial 
positive energy (/iMt)inj is extracted from the plasma and is deposited in the magnetic field to be carried 
outwards. Figures 10b and lib show the energy conversion between the fluid part and electromagnetic part 
in the Schwarzschild geometry, where Xfluid is always positive. Though the poloidal flow solution in the black 
hole magnetosphere contains two Alfven radii r = and r = r^"*, an accretion across both Alfven radii 
is possible when the injection point is located between the outer Alfven point and the outer light surface. 
One of them corresponds to the Alfven point for the considered flow, where the requirement of — is 
satisfied, and does not change its sign; such a point is A'" for case (i) and A°"* for case (ii). However, 
Xc-m changes its sign at the other Alfven radius, which is not the Alfven point for the considering SAF- 
solution because ^ there; such a radius is r = r™* for case (i) and r = r^* for case (ii). Outside this 
latter Alfven radius, B'^ /Bp < and magnetic energy streams outward {nu^\E\Xem > 0) to the injection 
point, while inside this point B'^/Bp > and magnetic energy streams inward. The fluid part of energy 
flux always streams inward ( r7,?i''|_B|A'[iuid < ), and it converts to electromagnetic energy flux as the fiow 
falls inward. If the ideal MHD plasma streams near the outer light surface, has a very large negative 
value, that means a very large outgoing magnetic energy flux in the flow. The magnetic field line is tightly 
wound up {—B'^/Bp ^1). To conserve the total energy flux along the accretion flow, the ingoing positive 
fluid energy flux should be also very large. When the poloidal motion of the plasma is very slow near the 
injection point, a large fluid energy flux must be due to the kinetic energy of the toroidal motion; so we can 
say that the origin of the large outward electromagnetic energy flux is the toroidal plasma motion near the 
plasma source. 

Figure 12 shows a negative energy accretion solution {ClpL > 1). We see that the Alfven point locates 
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inside the ergospherc (sec Paper I). The outgoing electromagnetic energy flux is always greater than the 
ingoing fluid energy flux {Xem < — 1 and Xfluid = — 1 — Xem > 0). The magnetic field lines are trailed 
{B'^/Bp < 0) everywhere due to the black hole rotation. 

For accretion with QpL < 1 and < < cou, the electromagnetic energy flux also streams outward 
everywhere, but at least near the event horizon the ingoing fluid energy flux dominates [i.e., — 1 < (Xem)H < 
0]. For magnetically dominated accretion, we see that X^m — —1. We should mention that the Poynting flux 
passing through the event horizon is not modified by the plasma inertia effect, explicitly. This is because the 
toroidal magnetic field at the event horizon becomes B^u = \/ (— (700/S)h (wh — ^F)ide'i')H for any ideal 
MHD accretion flows, which is the same expression as that of the force-free case (Znajek 1977). 

We should note that, in general, Cs Cf for an SAF-solution. So, when we try to estimate the locations 
of the fast and slow magnetosonic points for a given ?}, we can only obtain possible ranges of magnetosonic 
points by using the rcY-fj diagram. For example, if wc know the values of rp and Cf, we can find a possible 
region of rs for an acceptable Cs range. To determine the location of the slow magnetosonic point explicitly, 
we need to obtain the Cs value by solving the poloidal equation, which is a polynomial with high degree. 



4.3. Fluid-Dominated Flows 

The solution for case (ii) becomes hydrodynamical accretion in the hydro-dominated limit. In the limit 
of weak magnetic field (e = B^/ (SttP) <C 1), for a flow with Up{x) ~ 0{e°), we see that ~ 0(e~^); then, 
we obtain B^/Bp 0(e°), {E/n) = ut + 0{e), and {L/fj) = —u^ + 0{e). The poloidal equation (10) then 
becomes 



{l + ul) = -k(^^^ +0{e) 



(34) 



The function defined by equation (12) is reduced to A; = — gr** -|- 2g''^£ — g'f"f'£^, where £ = —u^/ut is the 
speciflc angular momentum as measured at infinity. The numerator (20) and denominator (21) become 



^^ = {M^Y(i + ui) 



Csw'5-^(l + 0(lnfc)' 



(35) 



V = {M'Y{Ct^-ul), (36) 

where the function <S = S{r,6) is related to the configuration of a stream line; for example, for a radial 

stream line, S{r) = — (2r^ — Smr + a^)/(rA). Here, we should remember that, in the ideal MHD case, a 
stream line coincides with a magnetic field line. So, even for a weak magnetic field, the factor iS should be 
related to the magnetic field lines, and in fact it must be replaced by 5 = (In Bp)' for MHD flows. Further, 
we should say that, in the weak-magnetic field limit, Q.f represents the angular frequency of the stream 
line for a hydrodynamical flow, and it would be determined as the angular velocity of the injection point. 
Thus, the above expressions are reduced to the relativistic hydrodynamical flow equations formulated by Lu 
(1986). In equations (24), (25) and (26), we can also check that the fast magnetosonic wave speed equals 
the sound wave speed, while the Alfven wave speed and the slow magnetosonic wave speed become zero. 

To the contrary, for ~ ^^aw ^ which corresponds to case (i), the hydrodynamical expression for 

the poloidal equation can not be obtained by a flow solution. At the Alfven point, M]^ = a. a ^ 0(e°) and 
u\ ~ 0(e). At the fast magnetosonic point, Up ~ 0(e). So, the velocity of trans- Alfvenic MHD accretion 
would also be the order of £^/^. 
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5. Concluding Remarks 

We have considered stationary and axisymmetric hot ideal MHD accretion along a flux-tube connected 
from a plasma source to the event horizon. To argue the details of the boundary conditions at the plasma 
source would take us beyond the scope of this paper. Therefore, we have surveyed the dependence of the 
trans- fast MHD flows on a wide range of source parameters. We have shown that, when the forbidden region 
is type lA or IIIA, there exist two physically different accretion regimes: (i) magneto-like MHD accretion and 
(ii) hydro-like MHD accretion. The magneto-like MHD accretion would be injected from near the separation 
point and passes through the inner Alfven point with a smaller 17. On the other hand, the hydro-like MHD 
accretion with a sufficiently large fj would be injected from between the outer Alfven point and the outer 
light surface and passes through the outer Alfven point. Hydro- like accretion may also be initially super- 
slow magnctosonic or supcr-Alfvcnic. A hot ideal MHD plasma with larger fj cannot accrete stationary onto 
the black hole after passing through the inner Alfven point. Then, if the value of fj increases with a secular 
timescale, the magneto- like MHD accretion solution should transit to the hydro-like MHD accretion solution; 
the inverse process would be also possible. The criterion for distinguishing between the two regimes is based 
on the locations of both the Alfven point and the fast magnetosonic point, which change discontinuously 
during the transition. 

For the magneto-like MHD accretion, we have found that the location of the X-type fast magnetosonic 
point is not unique for fixed intrinsic parameters of the accreting plasma. For example, in Figures 5a, 7b 
and 6, in the range of ru < rp < r^*, & fj = constant line crosses a solid curve of Cf < 0.1 at three points; 
the first and third fast magnetosonic points are X-type critical points, while the second one is 0-typc. So, 
we can expect two modes of magneto-like MHD accretion solutions; that is, one passes through the "inner" 
inner-fast magnetosonic point and the other passes through the "outer" inner-fast magnetosonic point. The 
number of these inner-fast magnetosonic points depends effectively on a, S, Qp, L, (^f and 77. There is a 
tendency for larger flp ~ luh to generate a multiple inner- fast magnetosonic solution. Compared with the 
radial field geometry, converging field geometries for accreting flows {S > 0) also generate such multiple inner- 
fast magnetosonic solutions. Transition between these two modes is also discontinuous. The ^-dependence 
of SAF-accrction solutions on the flow velocity u''{r) and the electromagnetic energy Xcmir) is very weak, 
as long as the location of the fast magnetosonic point does not jump, by changing the 6 value, to another 
branch of the multi- inner fast magnetosonic points. The radial terminal velocity at the event horizon u^, 
rather than the radial one, slightly decreases (increases) for ^ > (J < 0). For example, we have checked 
this for S = 0.4, 0.0 and —0.4 cases. For such accretion solutions, we also find that for (5 > ((5 < 0) the total 
energy flux per magnetic tube fjE and r) increase (decrease), while the total energy E decreases (increases). 

Throughout this paper, we have only discussed the ideal MHD flow cases. However, non-ideal MHD flow 
solutions near the event horizon are also presented (Punsly (1990, 2001)). Punsly (1990) discussed an ingoing 
magnetic flow solution along magnetic field lines that thread the ergosphere and the equatorial plane (and 
therefore not the event horizon). This solution corresponds to our sub-Alfvenic ingoing solution approaching 
the inner light surface with zero poloidal velocity (see, e.g.. Fig. 10a); along this solution B'^ = at r = r^* 
(see also, Fig. 10b), where Up ^ m^w (^o* Alfven point A), while for the SAF-solution B'^ 7^ at the 
Alfven point A. When we consider an accretion solution onto a black hole, it seems that there is no sub- 
Alfvcn accretion solution consistent with the ideal MHD approximations due to the existence of the forbidden 
region. However, of course, for such a set of ingoing flow parameters the ideal MHD approximation must 
be rejected, and then the non-ideal MHD ingoing flow should exist. This is because near the light surface 
due to the plasma inertia effects large radiation losses are expected, and large radiation losses equate to a 
dissipative plasma and a breakdown of the ideal MHD approximation. Then, a non-ideal MHD accretion 
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flow solution, which does not pass through the fast-magnctosonic critical point F discussed in the previous 
section, also exists in the region downstream of the light surface because of the inward attraction of black 
hole gravity. This requires that dissipative effects be incorporated into the physical description of the inward 
extension of the ideal MHD wind inside the light surface. Such a plasma propagates inside the inner light 
surface and enters the forbidden region of r < r™ with relatively slow velocity (as compared with the Alfven 
wave speed except for the area close to the horizon) by crossing or reconnecting the magnetic field lines, 
where the physical meaning of the concept of the forbidden region for the ideal MHD flows would be lost. 
Note that the non-ideal MHD accretion flow must also become super-Alfven and super-fast magnetosonic 
just inside the inner light surface (see Chapter 9 of Punsly 2001). 

For a "super-critical" accretion flow oi E > E-p and w^w < '^p < '^fmi which reaches the event 
horizon without passing through the fast magnetosonic point F and is unphysical under the ideal MHD 
approximation, some kinds of dissipative effects near the event horizon would also make the super-critical 
accretion onto the black hole possible. One may also expect a "sub-critical" accretion flow of E < Ep, which 
is also not a solution for the inflow into the black hole under the ideal MHD approximation. However, the 
breakdown of the ideal MHD would change the nature of the critical points (see, e.g., Chakrabarti (1990) for 
the hydro dynamical flow case); that is, instead of the X-type and 0-type critical points, the "nodal- type" 
and "spiral-type" critical points appear on the flow solutions in the (r, ■u'") -plane. Thus, when the initial 
condition at the plasma source does not match the critical condition at the fast magnetosonic point, the 
non-ideal MHD inflow into the black hole must be realized. The detailed structure of the non-ideal MHD 
flows around the critical point is a very important topic for black hole accretion, but further discussion is 
out of the scope of this paper. 

When the angular velocity of the magnetosphere is in the range of < ^Ip < wh, only one Alfven 
point appears on the r-u^ plane. To accrete onto the black hole, an ideal MHD flow must pass through the 
Alfven point classified by type HA*", IIA°"*, HB*" or IIB"' in Paper I. Determining whether the accretion is 
magneto-like or hydro-like may not be clear in type HA case. However, we can expect that the Alfven point 
A'" is responsible for the magneto-like accretion solution and the A"*** point is responsible for the hydro-like 
accretion solution. It seems that for the Alfven point of type HB""* there are no accretion solutions consistent 
with the ideal MHD approximation. Of course, for such a set of ingoing flow parameters the ideal MHD 
approximation must be abandoned. When the SAF-accretion solution is invalid (e.g., the type HE""* case), 
the non-ideal MHD accretion becomes important and would dynamically effect the magnetic fleld structure. 

When a region of k(r) > exists along a magnetic field line for larger L^lp (< 1), the forbidden region 
is type B and the ideal MHD accretion is only allowed after passing the inner Alfven point. The hydro-like 
MHD accretion does not arise because of the sufficiently strong centrifugal barrier, while the magneto-like 
MHD accretion is available because of the effective angular momentum transport from the fluid-part of 
total angular momentum to the magnetic-part (Hirotani ct al. 1992). When the hot effects dominate in the 
plasma, this magneto-like MHD accretion may be also forbidden due to the disappearance of the inner-fast 
magnetosonic point. In this case, however, we can also expect the non-ideal MHD accretion to fall into the 
black hole. 

If a shock front is generated after passing the fast magnetosonic point, the post-shock flow with increased 
entropy must pass another fast magnetosonic point again on the way to the event horizon. To construct 
such a shock formation model, the existence of multiple fast magnetosonic points is required in the accretion 
solution. We can expect two types of discontinuous transitions. One is the transition from the hydro-like 
solution to the magneto-like solution at somewhere between the middle-fast and inner-fast magnetosonic 
points. The other is the transition in the hydro-like or magneto-like solution. For the magneto-like solution. 
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we can see the possibility of transition from the rcT-fj diagrams. That is, the accreting matter passes 
through the first inner-fast magnetosonic point located just inside the inner-Alfven point, and then after 
the shock it passes through the third (X-type) inner-fast magnetosonic point. Here, the first and third 
inner-fast magnetosonic points are located on different Cr-curves in the ra-fj diagram because of the entropy 
generation. The insights gained in the course of our analysis should be of use in further investigations of 
shocked accretion solutions. 

Although we have discussed accretion flows onto a black hole, our treatment can be applied to outgoing 
flows (i.e., winds and jets). To do so, we can plot Tcr-V diagrams in the range of < rcr < oo (for the fast 
magnetosonic point, < rp < oo; for the slow magnetosonic point, r™ < rs < r^"*). So, we will find the 
possible locations of the magnetosonic point for outflows. When the forbidden region is type lA or IIIA, we 
flnd two regimes of SAF-solutions for outflows; that is, inj S™"* — * A°"* — » F°"* — > oo and inj — » S™ —^ A" 
_^ pmid _^ where the slow magnetosonic point S™ is located between the inner light surface and the inner 
Alfven point, and the fast magnetosonic point F°"* is located outside the outer Alfven point. The former 
SAP-solution is also available for the type B forbidden region, and the latter SAF-solution remains in the 
hydrodynamical limit. From the Tcr- ry diagrams, for the outflow, the location of the fast magnetosonic points 
with the same fj value strongly depends on 6. A similar result has been discussed in Takahashi & Shibata 
(1998) as a pulsar wind model without the gravitational effects. Here, we should note that to blow away to 
infinity, the outgoing flows also must pass through the slow magnetosonic point and the Alfven point, but at 
the asymptotic region both the super-fast and sub-fast magnetosonic outflows are available. Whether or not 
the outflow passes through the fast magnetosonic point depends on the field aligned parameters of flows and 
the geometry of fleld lines. So, for the outflow, the fast-magnetosonic surface does not need to cover all solid 
angles. We can expect that the trans-fast MHD outflow is realized at least in some part of the magnetic field 
lines, and the distribution of the fast-magnetosonic surface would play a very important part in explaining 
the generation and coUimation of a highly accelerated jet or wind. 

Wc arc grateful to Akira Tomimatsu, Sachiko Tsuruta and Vasily S. Bcskin for useful discussions and 
to an anonymous referee for constitutive criticism which helped us to improve the paper. 
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Fig. 1. — MHD accretion onto a black hole. The flow injected from a plasma source (e.g., the disk surface or 
the corona) passes through the slow magnetosonic point, the Alfven point and the fast magnetosonic point, 
in order, and then goes across the event horizon. 




Fig. 2. — Schematic pictm-es of V{r,u'^) — curves for (a) strong magnetic field (Cg^)A < u\ and (b) weak 
magnetic field {C^^)a > u\. The thick-solid-curves correspond to the — u^yi curves; the thick-dash-dot- 
curves correspond to the = u^w curves; and the thick-dotted-curves correspond to the = Ugj^j curves. 
The shaded regions are forbidden regions of flows (type lA or IIIA), and separate the super- Alfvenic region 
labeled by "super" from the sub- Alfvenic region labeled by "sub". The hatched region is inside the black 
hole. 
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Fig. 4. — (a) Relations between rcr and i) and (b) relations between E/ jic and for a radial flow with 5 = 
in a Schwarzschild black hole magnctosphcrc with several Q,- values (0.00, 0.10, 0.20, 0.30). The angular 
velocity of the magnetic field is 0_f = O.Sf^max, where milmax = 0.192. The location of the Alfven point 
is specified by xa = 0.8, which corresponds to QpL = 0.7287. The solid-curves indicate fj = fj{rF) with 
7'ci = ^'F- The dashed-curves indicate fj = f){rs) with rcr = rs- The shape of forbidden regions corresponds 
to type lA. 
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Fig. 5. — Relations between 17 and for a radial flow with 5 = in a Schwarzschild black hole magnetosphere 
with several values (0.00, 0.10, 0.20, 0.30). The angular velocities of the magnetic field are (a) ftp = 
0.9<linax and (b) Of = O.SOmax- The location of the Alfven point is specified by xa = 0.8, which corresponds 
to (a) i^pL = 0.8841 and (b) QpL = 0.2687. The solid-curves indicate fj = fjirp), and the dashed-curves 
indicate fj = fj{rs)- The shape of the forbidden regions corresponds to type lA. 
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Fig. 6. — Relations between f] and r„ for a radial flow with 5 = and several Ccr values (0.00, 0.10, 0.20, 0.30) 
in a Schwarzschild black hole magnetosphere. The angular velocity of the magnetic field is Of = O.SOmax- 
The location of the Alfven point is specified by xa = 0.5, which corresponds to VtpL = 0.6434. The solid- 
curves indicate r) = fjirv), and the dashed-curves indicate fj = ry(rs). The shape of the forbidden regions 
corresponds to type lA. 
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Fig. 7. — Relations between r) and rcr for a radial flow with (5 = in a Kerr black hole magnetosphere; 
(a) a = 0.5m and (b) a = —0.5m with several Q„ values (0.00, 0.10, 0.20, 0.30). The angular velocity of 
the magnetic field is Q.f = O-Silmax, where (a) mfimax = 0.244 and (b) mfimax = 0.163. The location of 
the Alfven point is specified by .ta = 0.8, which corresponds to (a) Q.pL = 0.6749 and (b) VlpL = 0.7490. 
The solid-curves indicate i) = f}{r-p), and the dashed-curves indicate r) = fi{rs)- The shape of the forbidden 
regions corresponds to type lA. 
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Fig. 8. — Relations between fj and rcr for a radial flow with 5 = and several Ccr values (0.00, 0.10, 0.20, 0.30) 
in a Kerr black hole magnetosphere (a = 0.8m). The angular velocity of the magnetic field is Q.f = O.SOmax, 

where mfimax = 0.309. The location of the Alfven point is specified by .ta = 0.8, which corresponds to 
Q,fL = —0.0117. The solid-curves indicate i) = fi{rF), and the dashed-curves indicate 77 = r?(rs). The shape 
of the forbidden regions corresponds to type IIA. 




Fig. 9. — The ^-dependence {5 = —0.4,0.0,0.1,0.4) of the relations between fj and rcr, where a = 0.0, 

Ccr = 0.1, r^F = O.Sr^max and XA = 0.8 (sec also Fig. 4a). The solid-curves indicate fj = ??(rF), and the 
dashed-curves indicate fj = r?(rs). The shape of the forbidden regions corresponds to type lA. 
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Fig. 10. — An example of the trans-fast MHD accretion solution (thick curves); (a) radial 4-velocity and (b) 
the ratio of the electromagnetic energy to the total energy. This solution passes through the inner Alfven 
point and the inner fast magnetosonic point (S ^ A ^ F ^ H), where a = 0.0, flp = O.SJlmax, S = 0.0, 
r = 4/3 and (p = 0.2. The location of the Alfven point is specified by xa = 0.8, which corresponds to 
QpL = 0.7287. The location of the fast magnetosonic point is rp = 2.10935m, which gives E-p/ ji^ = 3.6643 
and 77 = 0.0337. The requirement that i) < ?7max is satisfied. 
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Fig. 11. — An example of the trans-fast MHD accretion solution (thick curves) passing through the outer 

Alfvcn point and the middle-fast magnctosonic point; (a) radial 4- velocity and (b) the ratio of the electro- 
magnetic energy to the total energy. The location of the fast magnetosonic point is rp = 3.00490m, which 
gives -Bf/Mc = 9.8666 and fj = 0.0378. The other parameters are the same as in Fig. 10. 



-27- 




Fig. 12. — An example of a negative energy trans-fast MHD accretion solution (thick curves); with (a) radial 
4-velocity and (b) the ratio of the electromagnetic energy to the total energy. The flow parameters are 



0.8m, Qp = O.ea, 



0.0, r 



which gives E-p/ ji^ = —0.15855 and i) 



= 4/3, Cf = 0.1, XA = 0.95 {Q.fL = 5.1710) and rp = 1.647934m, 
0.0297. 



